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Abstract 

First-principles calculations within the DFT are employed to investigate the relationship between 
the electronic structure and the unexpected features of the hexagonal cell parameters of YbAuIn 
and Yb3AuGe2ln3. Calculations indicate that YbAuIn is an intermediate valent system with one 
Yb 4 / state pinned to the Fermi level, while Yb3AuGe2ln3 is closer to integer valency with all 
Yb 4 / states occupied. Structural relaxations performed on LaAuIn and LuAuIn analogs reveal 
that expansion of the c-parameter in Yb3AuGe2ln3 is attributable to larger size of the divalent Yb 
compared with intermediate valent Yb. 
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I. INTRODUCTION 


Among strongly correlated materials, Yb-containing intermetallics are particularly fas¬ 
cinating, because they display a variety of intriguing physical properties arising from the 
presence of the localized 4/ states 1 . Yb exhibits two valence states: magnetic, trivalent 
(Yb 3+ ) or nonmagnetic, divalent (Yb 2+ ) configurations. Depending on the electronic and 
crystallographic environment, as well as on external conditions of temperature and pressure, 
the two possible valence configurations can be energetically degenerate, giving rise to mixed 
and intermediate valency, magnetic transitions induced by pressure or alloying, anomalous 
low temperature thermodynamic and transport propertied- 23 . Furthermore, because of the 
hybridization between the localized 4/ states and delocalized s and p states many Yb-based 
intermetallics exhibit heavy fermion behavior, Kondo coupling, Ruderman-Kittel-Kasuya- 
Yosida interaction 1 2 4 11 . Families of isostructural Yb-containing systems display a wealth 
of ground state properties and different low energy excitations as a function of composition. 
Far from being exhaustive, examples of such systems include YbMCuj (M=Ag, Au, Pd, In, 
Cu, Cd, Mg, Tl, Znj™, YbTX (^transition metal, X=Sn, an d Yb 2 T 2 In (T=Cu, 

Pd, Au)P. 

The present theoretical investigation was motivated by a recent experimental study on the 
crystal structure, electronic and magnetic properties of the newly synthesized Yb 3 AuGe 2 In 3 
and its ternary analog YbAuIn 24 . As shown in fig. [lj the new compound Yb 3 AuGe 2 In 3 has 
a hexagonal crystal structure (P62m) and it is obtained by replacing two Au atoms in the 
structure of YbAuIn 25 ' with two Ge atoms. 

The crystallographic analysis of the two Yb-containing systems reveals a peculiar fea¬ 
ture 24 : when the two Au atoms in YbAuIn are replaced by Ge, the hexagonal unit cell 
exhibits a contraction along the a and b axes and an expansion along the c axis. On the 
one hand, the fact that the a- and ^-parameters of YbAuIn ( a=b= 7.712 A) are larger than 
those of Yb 3 AuGe 2 In 3 (a = b = 7.315 A) makes sense because the Au atom is larger than 
Ge, but on the other hand, the observation that the c-parameter of YbAuIn (c = 4.029 A) 
is 10% shorter than that of Yb 3 AuGe 2 In 3 (c = 4.421 A) cannot be explained by the size 
difference between Au and Ge. The present work addresses the puzzle of this unusual struc¬ 
tural change using electronic structure calculations on YbAuIn and Yb 3 AuGe 2 In 3 , focusing 
on the role played by the Yb 4/-electrons in the electronic and structural properties of the 
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FIG. 1. (Color online)Hexagonal unit cell of Yb 3 AuGe 2 ln 3 . The structure cen be regarded as 
alternating monoatomic layers of In 3 Ge 2 (or In 3 Au 2 in the case of YbAuIn) and YbsAu. 

two systems. 

II. METHOD OF CALCULATION 

The calculations have been carried out using the highly accurate full potential linearized 
augmented plane wave method (FP-LAPW) 26 within density functional theory (DFT) 2 ' 28 
as implemented in the WIEN2k program 29 . The exchange and correlation potential was es¬ 
timated by the local spin density approximation (LSDA) 30 , necessary for magnetic systems. 
To overcome the inability of DFT to model localized states, on-site Coulomb interaction was 
added to the Yb 4/ states using the “around-mean field” version (AMF) of the LSDA+ U for- 
malismf 31 ! According to the literature physically reasonable values for the screened Coulomb 
repulsion (Hubbard U parameter) within the Yb 4/ states are in the range of 6-7 eV. To 
check the dependence of the results on the particular value of the U parameter, the cal¬ 
culations have been performed for U = 5.44, 6.75, and 8.16 eV (corresponding to 0.4, 0.5, 
and 0.6 Ry). Based on energy convergence tests, the calculations were carried out using the 
following setup: the muffin-tin radii values (in atomic units; 1 au = 0.529 A) were chosen as 
2.5 for all the atoms, the cutoff parameter was RK max = 9.5, and the energy of separation 
between valence and core states was —6.0 Ry. Self-consistent iterations were performed with 
35 k-points in the irreducible part of the Brillouin zone (IBZ), and convergence was assumed 
when both the energy and the charge difference between two consecutive iterations were less 
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FIG. 2. (Color online) (a) Band structure of Yb 3 Au 3 lii 3 with the Au s orbital character emphasized. 
The partially occupied 4/ level lies at E^r, while the occupied levels are spread around -1.0 and 
-2.4 eV. (b) Total (shaded area) and partial density of states of YbAuIn associated with the Yb 4/ 
orbitals, calculated for different values of Hubbard U. The partially occupied / state is pinned at 
the Fermi level. 


than 0.0001 Ry and O.OOle respectively. Scalar relativistic corrections were added and spin- 
orbit coupling (SOC) was incorporated using second variational procedure 32 . With SOC 
included, the energy range in which eigenvalues were searched, was set from —9.0 to 7.0 Ry. 
The calculations were carried out using the experimentally determined lattice parameters. 


III. RESULTS AND DISCUSSIONS 

The calculated band structure of YbAuIn is illustrated in fig. |2](a) and reveals the interme¬ 
diate valent nature of YbAuIn, in agreement with magnetic susceptibility measurements 3 24 . 
The flat bands located around 1.0 and 2.4 eV below the Fermi level (E F ) correspond to the 
occupied 4/ states of the Yb ion. These states are split into 4/ 5 / 2 and 4 / 7/2 manifolds under 
spin-orbit interaction of approximately 1.3 eV. Each of the two spin-orbit complexes are fur¬ 
ther split due to the anisotropy of the Coulomb interaction within the 4/ shell and possibly 
due to crystal field effect! 33 . The additional flat band located at the E F corresponds to a 
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partially occupied 4/ level. It is remarkable that the pinning of the 4/ state at the Ep occurs 
even though, through manipulation of the density matrices^, the calculations were initialized 
with the Yb 4/ shell fully occupied. Apparently, during the process of self-consistent relax¬ 
ation one Yb 4/ band is split off from the rest of the 4/ complex, becomes partially occupied 
and is pinned to the Ep. This type of electronic structure is characteristic to intermediate 
valent heavy fermion compounds and it has been observed in other Yb-containing systems, 
such as YbBiPtpS and Yb 4 As^51. 

Figure |2^b) illustrates the spin polarized density of states (DOS) of YbAuIn, projected 
on the Yb 4/ orbitals for different values of U. As the U -parameter increases the partially 
occupied 4/ level remains pinned at the Ep while the occupied levels shift down on the 
energy scale. The occupied and the partially occupied 4/-states are well separated, indicating 
intermediate valent character for Yb in YbAuIn. To evaluate the valency of Yb the traces 
of the density matrices have been computed. The calculated 4/ occupations as a function 
of the [/-parameter indicate valences of 2.56, 2.60, and 2.66 for Yb for U = 5.44, 6.75, and 
8.16 eV respectively, somewhat larger that the experimental valence of 2.1« However, 
if we consider the fact that in the LAPW method the partial charges are considered only 
within the muffin-tin radii, and there is charge left in the interstitial region, the calculated 
valences are in reasonably good agreement with the experiment. 

In the case of YbsAuGe 2 ln 3 the calculations indicate a quite different picture for the 
electronic structure and Yb valency. The band structure of YbsAuGe 2 ln 3 , illustrated in 
fig- 3 a) displays two flat bands located at approximately 0.2 and 1.5 eV below Ep. These 
are the Yb 4/ bands, split into 4/ 5 / 2 and 4 f 7 / 2 manifolds, separated by the 1.3 eV spin orbit 
interaction. This can also be observed in fig. [3](b) , which illustrates the DOS associated 
with the Yb 4/ levels, calculated for different [/-parameters. It is apparent that the Yb 4/ 
states are occupied, located below Ep and their position relative to Ep does not change for 
different values of the Coulomb repulsion U. 

Because YbAuIn and Yb 3 AuGe 2 ln 3 are isostructural and the latter is obtained by re¬ 
placing two Au atoms in YbAuIn (or Yb 3 Au 3 ln 3 ) by two Ge atoms, the differences in the 
electronic structure must be related to the presence of Ge. The physics of the Yb 4/ occu¬ 
pation can be understood by analyzing the positions of the 4/ levels relative to the Au 6s 
bands in YbAuIn and Ge 4 p in Yb 3 AuGe 2 hi 3 , as illustrated with fat band representations 
in figs, ga) and[3](a). In YbAuIn the Au 6s states lie mostly below the Yb 4/-level. When 
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FIG. 3. (Color online) (a) Band structure of Yb 3 Au 3 lii 3 with the Au s orbital character emphasized. 
The partially occupied 4/ level lies at Eyr, while the occupied levels are spread around -1.0 and 
-2.4 eV. (b) Total (shaded area) and partial density of states of YbAuIn associated with the Yb 4/ 
orbitals, calculated for different values of Hubbard U. The partially occupied / state is pinned at 
the Fermi level. 


Au (5d 10 6s 1 ) is replaced by Ge (3d 10 4s 2 4p 2 ) the number of valence electrons in the system 
increases, therefore the Fermi level is shifted up in energy. Because the Ge 3 p band is not 
completely filled, the bands near E^(both below and above Ep) display significant Ge p 
character. Therefore, the upper Yb 4/-level, which was initially pinned at the Ep, now 
lies below the top of Ge 3 p band, becomes fully occupied just below the E F . Therefore, 
according to electronic structure calculations that include magnetic and relativistic effects 
as well as intra-atomic correlation, there is an essential difference between the two systems: 
in the YbAuIn the Yb ion has intermediate valent configuration, with one partially hlled 4/ 
level and a valency of ~2.6, while in Yb 3 AuGe 2 lii 3 , Yb is calculated to be closer to divalent 
state (Yb 2+ ). 

The unusual behaviour of the lattice parameters, observed experimentally®, can be asso¬ 
ciated with the valence configuration of the Yb. Because Ge is smaller than Au, one would 
expect overall shorter lattice parameters forYb 3 AuGe 2 ln 3 compared to those of YbAuIn. 
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However, as shown above, in the presence of Ge, the valency of Yb becomes closer to Yb 2+ 
and because divalent Yb (Yb 2+ ) is larger than a mixed valent Yb (Yb 2 6+ ), the overall 
decrease in the lattice constants due to Ge is be counterbalanced by the increase in the 
size of Yb. If we consider the crystal structure of YbAuIn as being formed by alternating 
monoatomic sheets of I 113 AU 2 and YbsAu stacked along the c-axis, we have the following 
scenario: on the one hand the area of the I 113 AU 2 sheets shrinks when Au is replaced by Ge, 
but on the other hand, the separation between them increases, because the size of the Yb 
ion increases as the result of Yb 2 ' 6+ —> Yb 2+ valence transition. The net effect is a decrease 
in the a- and ^-parameters and an increase in the c-parameter of the hexagonal unit cell. 

To verify this assumption, full structural relaxation have been carried out on two hy¬ 
pothetical compounds, LaAuIn and LuAuIn where La and Lu (La > Lu) were chosen to 
simulate the divalent and intermediate valent Yb, respectively, and to avoid complications 
arising from the presence of localized / electrons. The calculations were performed using 
VASP package 36 37 and included relaxation of internal parameters (atomic positions) and 
volume optimization. The calculated lattice parameters for LuAuIn are a = 7.79 A, c = 
3.85 Aand for LaAuIn are a = 7.89 A, c = 4.32 A, that is, when the smaller Lu is replaced by 
the larger La the a-parameter increases only by ~1.3%, while the c-parameter increases by 
~12.2%. This indicates that the size of the ions that are located between the monoatomic 
sheets of In 3 Au 2 primarily influences the length of the c-parameter in the hexagonal unit 
cell, and has only a minor effect on the a- and ^-parameters. This supports the premise 
that the unusual structural relationship between YbAuIn and Yb 3 AuGe 2 hi 3 , is associated 
with the Ge-induced Yb 2 ' 6+ —> Yb 2+ the valence transition. 

In summary, using first-principles theoretical methods, the electronic structures of two 
isostructural Yb-intermetallics, YbAuIn and Yb 3 AuGe 2 ln 3 , have been investigated, focusing 
on the valence configuration of the Yb ion. YbAuIn is calculated as intermediate valence 
compound, with one partially filled Yb 4/ state, pinned to the Fermi level, which sug¬ 
gests that YbAuIn is a heavy fermion compound. Yb 3 AuGe 2 ln 3 is calculated to be closer 
to integer valency with a completely occupied Yb 4/ shell. Furthermore, the electronic 
structure calculations combined with structural optimizations performed on LaAuIn and 
LuAuIn analogs, indicate that the longer c-parameter in Yb 3 AuGe 2 ln 3 compared with that 
in YbAuIn is attributable to the fact that the divalent Yb ion in Yb 3 AuGe 2 ln 3 is larger than 
the intermediate valent Yb in YbAuIn. The change in the valence configuration is caused 
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by the presence of Ge. 
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